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Abstract 

> 

This note aims to subsume several apparently unrelated models under a common frame- 
work. Several examples of well-known quantum field theories are listed which are con- 
nected via stochastic quantization. We highlight the fact that the quantization method used 
to obtain the quantum crystal is a discrete analog of stochastic quantization. This model 
is of interest for string theory since the (classical) melting crystal corner is related to the 
topological A-model. We outline several ideas for interpreting the quantum crystal on the 
string theory side of the correspondence, exploring interpretations in the Wheeler-De Witt 
framework and in terms of a non-Lorentz invariant limit of topological M-theory. 
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1 Introduction 

In this note, we aim to connect several apparently unrelated theories by showing that they 
are related through a common quantization scheme, namely stochastic quantization [1J. It 
finds applications in field theory, the world of statistical and condensed matter physics, and 
via correspondences also in string theory. 

Models that can be regarded as resulting from stochastic quantization all have a common 
structure. Starting from a (classical) model in d dimensions one obtains a quantum model in 
d + 1 dimensions as a Markov process that converges at thermal equilibrium to the minima 
of the classical action. The partition function of the classical precursor model is expressed 
as the norm square of a wave function, which is the ground state of the theory in d + 1 
dimensions, 

Ztl = (^undl^ound) • (I'D 

After summarizing the preliminaries of stochastic quantization, we collect several exam- 
ples of theories which are related by stochastic quantization. The most basic and prototypical 
example is the one of zero dimensional field theory, whose stochastic quantization leads straight 
to supersymmetric quantum mechanics. In the case of our next example, the free boson, the 
naturally supersymmetric description of the stochastically quantized theory provides a true 
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advantage over its prior treatment in the literature Q. We show furthermore that the quan- 
tization scheme employed to obtain the quantum crystal in [3], which is the same as the one 
used for the quantum dimer model [4], is nothing else but a discrete analog of stochastic 
quantization. In fact, both schemes are based on an underlying Markov process, which is a 
stochastic process for which the probability for the system to be in a state at time t depends 
only on the state at t — 1. 

Even though stochastic quantization cannot be applied directly to the gauged wzw model, 
its partition function can be expressed as the norm square of a 3-dimensional wave func- 
tion, which is the ground state of a strong coupling limit of topologically massive gauge 
theory (Chern-Simons plus Yang Mills). This example thus belongs to the same class of 
models. The example of the quantum crystal is interesting from the point of view of string 
theory, since a correspondence between the classical crystal melting configurations and the 
topological string A-model has been found (5l IQ. The quantum crystal, giving rise to a 
seven dimensional theory, suggests a connection to topological M-theory. We outline some 
ideas concerning the interpretation of the quantum crystal on the string theory side of the 
correspondence, such as a Wheeler-De Witt interpretation of the eigenvalue equation for the 
ground state. Inspired by the crystal/ string correspondence, also the stochastic quantiza- 
tion of Kahler gravity [7] is sketched. Since Kahler gravity is the target space description of 
the topological A-model, its quantization should be related to the quantum crystal, and via 
Hitchin's functional also to topological M-theory (8|. Given that quantum theories such as 
those arising from the stochastic quantization of Kahler geometry are all in need of regular- 
ization, the discrete framework of the quantum crystal ultimately seems, if anything, most 
suited for practical purposes. 

The plan of this note is as follows. In Section |2j we summarize the technique of stochas- 
tic quantization in the Langevin and Fokker-Planck formulation and present its manifestly 
supersymmetric form, as first given in H9l. In Section 2.4 the discrete analog of the stochastic 



quantization scheme is discussed. In Section |3j several examples of theories that are con- 
nected via stochastic quantization are given. The prototypical example is the one of zero 
dimensional field theory, discussed in Section |3.1| The next example is the free boson, see 



Section[32| In Sec. 3.3 we discuss the relation of the gauged wzw model to a strong coupling 



limit of topologically massive gauge theory. In Sec. |3.4[ we briefly mention the quantum 



dimer model which is a precursor of the quantum crystal. The latter is treated in Sec. 3.5 



where also possible interpretations of the quantum system in string theory are outlined. The 



last example concerns the stochastic quantization of Kahler gravity, see Sec. 3.6 A brief 
summary and conclusion is given in Sec.|4] 
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2 Stochastic quantization revisited 



Stochastic quantization is a quantization method for Euclidean field theories introduced in 
1981 by Parisi and Wu [lj. It makes use of the fact that Euclidean Green's functions can be 
interpreted as correlation functions of a statistical system in equilibrium. A good overview 
over this quantization scheme including also subsequent work can be found in IfTUI . 

Starting from a d dimensional Euclidean field theory the field is coupled to a white Gaus- 
sian noise which forces it to a random movement on its manifold, a continuum analog of 
a random walk. This stochastic system evolves along a new, fictitious, Euclidean time t. 
In the limit t —*■ oo, where thermal equilibrium is reached, the d-dimensional correlation 
functions of the quantum field theory are recovered, thus the original d-dimensional field 
theory has been quantized. The equivalence of stochastic quantization and conventional 
path integral quantization has been shown. When considered purely as a way to quantize 
a d-dimensional field theory, the extra time dimension is but a computational device which 
allows us to compute the correlators of the d-dimensional quantum field theory after the 
id + 1) -dimensional system has settled into equilibrium. In this paper, we will study the 
(d + 1) -dimensional quantum theory in its own right. 

In the following, we will give a very brief introduction to the subject, referring the reader 
to the literature for the formal definitions of all the concepts used. We follow the notation 
used by GO). 

2.1 Langevin formulation 

A stochastic variable X, defined by a range of values x and a probability distribution P(x) over 
these values, is generally used to describe the random fluctuations of a heat reservoir back- 
ground. A stochastic process is a process depending on the time t and a stochastic variable. A 
Markov process is a stochastic process in which the conditional probability for the system to 
be in a certain state at time t only depends on its state at t — St and not on its earlier history. 
The simplest example of a Markov process process is Brownian motion. Take a particle of 
mass m moving in a liquid with friction coefficient a. Its motion obeys the Langevin equation 

m^-v = —ocv(t) + fj{t), (2.1) 
at 

where fj(t) is the stochastic force vector representing the collisions between the particle and 
the molecules of the fluid. fj(t) is a stochastic variable and can be assumed to have a Gaus- 
sian distribution. 

Stochastic quantization of a Euclidean field theory works as follows. We supplement the 
field (p(x) with an extra time dimension t (which must not be confused with the Euclidean 
time Xq). Then we demand that the time evolution of (p(x, t) obeys a stochastic differential 
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equation such as the Langevin equation ( p.!} , which allows the relaxation to equilibrium: 

d(p(x,t) 1 SS C \ 



dt 2 Sep 



+ r](x f t), (2.2) 



with S c i the Euclidean action. The correlations of r/, which is a white Gaussian noise, are 
given by 

(rj(x / t))=0 > {rj{x 1 ,t 1 )r ] {x 2 ,t t )) =25{t 1 -t 2 )S d {x 1 -x 2 ). (2.3) 



Equation \2.2\ has to be solved given an initial condition at t = to leading to an ^/-dependent 
solution (p,j (x, t). As a consequence, also (x, t) is now a stochastic variable. Its correlation 
functions are defined by 

The central points for stochastic quantization to make sense is that equilibrium is reached 
for t — > oo, and that 

lim (^ t ,(x 1 ,t)...(p r j(x k ,t)) = {<p(x 1 )...<p(x k )) , (2.5) 

i.e. that the equal time correlators for <p n tend to the corresponding quantum Green's func- 
tions. 

2.2 Fokker-Planck formulation 

The conditional probability distribution of the speed of a particle executing a Brownian mo- 
tion depends on the initial conditions P(v, to) and satisfies the Fokker-Planck (FP) equation: 



which can be derived from the Langevin equation ( |2.1| | in a standard way (after setting 
m = a = 1). For stochastic quantization, the FP formulation can be used as an alternative 
approach. 



Starting from (2.2 \, the evolution of the probability P[(p(x), t] in the extra time direction t 



is given by the FP equation 
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together with an initial condition. One easily verifies that for t — > oo, 

The FP equation is a so-called master equation (which means in particular that its kernel is 
not Hermitian). To bring it into Hamiltonian form we introduce 

ip[p,t] = P[(p,t]e Sd ^ n . (2.9) 

This function tp satisfies the equation 

M, *\ = - / d d xH[cp} t/#, t] , (2.10) 



a 

at' 



where 

«W = -|^ + "W. and UW = I(^) 2 -l^!. (2.11, 
The zero energy ground state for this Hamiltonian is 

tp [$] = e -SdW/2_ (2.12) 

The corresponding Lagrangian can be obtained as usual by a Legendre transform: 

£[<P,<P} = 1<P 2 + U[<p}. (2.13) 

An alternative and more general way to find the quantum action consists in considering 
the partition function 

Z = = J Vrje-i^ 2( - x ^ dxit . (2.14) 
Changing variables from rj to (p, this becomes 

Z = jv<p det (|) r^(^ + W) V (2.15) 

The evaluation of the Jacobian requires some care. In particular, one has to choose a direc- 
tion for the time evolution, make use of the Stratanovich stochastic calculus (see flTTJ), and 
drop an infinite constant. All these issues are neatly resolved in the supersymmetric formu- 
lation presented in next section. By choosing propagation in the positive time direction, one 
formally recovers 



det ( j^j = e u s <p 2 . (2.16) 



Plugging the Jacobian back into (2.151 and neglecting a total derivative term (pSS = S, one 
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finds 



where 



Z = [ Dcpe- f c ^ dxdt , (2.17) 



2 r 8 \ Sep J 4 Sep 2 ' 



as in Eq. ( 2.13| . 



2.3 Supersymmetric formulation 

It was shown in [9] that any theory obtained by stochastic quantization admits supersym- 
metry with respect to the fictitious time direction t. In fact, already the classical Brownian 
motion in one dimension corresponds to supersymmetric quantum mechanics. This super- 
symmetry is linked to the existence of the so-called Nicolai map. One of the useful con- 
sequences of this supersymmetry which is already encoded in the Langevin equation, are 
Ward identities for the Green's functions. 

The fermionic superpartners of <p can be introduced as a calculational device to expand 
the Jacobian determinant ( 2.16| as a fermionic path integral: 



The FP Lagrangian becomes 

CM^^^ + l^-^l + ^y (2.20, 

In the scalar case, this can be rephrased in terms of superspace language. We introduce the 
superfield 

O = cp + Sip + xpe + 66F, (2.21) 

where F is an auxiliary field needed to achieve closure of the supersymmetry algebra. The 
super-covariant derivatives are given by 

With this, the action becomes 

C FP [O] = S d [O] + ODDO, (2.23) 

where S c \ [O] is the integral over the classical Lagrangian in terms of the superfield O. Note 
in particular that S c \ [O] now takes the role of a superpotential which is a function of the su- 
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perfield O. A word should be spent on the interpretation of the supersymmetry in this type 
of problems. When passing from the Langevin to the FP description, there are two possible 
choices for the direction of the time propagation. They correspond to two Hamiltonians, 
Or and H + . H~ is the one we used; TL + has by construction only strictly positive eigen- 
values and hence does not contribute to the t — » oo dynamics where the classical system is 
recovered. The partition function can be expressed as the difference Z = Z + — Z using the 
corresponding Lagrangians. For some functions F[(p] of the fields, the expectation values on 
the two partition functions coincide and we find identities of the type 



(F)z = <F) Z+ - (F) z - = 0, 
which can be interpreted as Ward identities in the supersymmetric formalism. 



(2.24) 



Consider now the special case of the action for a free field <p which can be written in the 
form 

SM = (<p\L<p) , (2.25) 

where <p G , with Jf 7 being a Hilbert space with scalar product (• |-) and L a self-adjoint 
operator. The Langevin equation for <p t reads 



<pt = -L(f> t + t] t , 

where t]t G ^ is a white Gaussian noise. The corresponding partition function is 

"1 



(2.26) 



Z = J T>t]t exp 
Changing the integration variable, 



dt (t] t \t]t 



(2.27) 



z = /^' del (t) e "p[i/* 



(2.28) 



We introduce the Grassmann fields xp,xp G J$? to express the determinant in a fermionic 
representation, 

(2.29) 



(2.30) 



Now the final action takes the form 



s^uMt4t\ = I \\<P + l<P\\ 2 - ($\ip + Ly) ■ 



This same action can be written in a manifestly supersymmetric form as 

S q [0] = J d0d6S susy = dM6 (DO|DO) + (0|LO) 



(2.31) 
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2.4 Discrete analog 

It is possible to quantize discrete stochastic models in a manner analogous to stochastic 
quantization as discussed above. The dynamics of any discrete stochastic model is described 
by the master equation 

A Pa (t) = £ (W a pPp(t) - Wp a P a (t)) , (2.32) 
at p 

where P a (t) is the probability to be in configuration a. at time t, and W a p is the transition rate 
to state a if the system is in state f>. We choose the following transition rates: 

W^ = C^-*(*M-*<«) /2 , (2.33) 

where g is a coupling constant, is the classical Hamiltonian evaluated on the configu- 
ration a and C a p is the adjacency matrix of the state graph. One can verify that the (unique) 
stationary distribution for P a is 

pf = 1 -M«) t (2.34) 

where 

Z = £g-««W. (2.35) 



This is equivalent to saying that the system is in a Boltzmann distribution with classical 
energy H. One can easily verify not only t 
also satisfies the detailed balance condition 



energy H. One can easily verify not only that P^ describes an equilibrium state, but that it 



WpA 0) =Wa?PP- (2.36) 

Detailed balance implies the Markov property for a stochastic process. 
It is customary to define the exit rate from state a. as 

W aa = -E W P« = -L C afi e-^-^V . (2.37) 
In this way, the evolution can be described by the vector equation 

^P(0=WP(0, (2.38) 
where W is the matrix with entries W a p. The stationarity condition becomes 

WP(°> = . (2.39) 
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Instead of using the matrix W, one can define the symmetrized version with entries 

W afi = (pi 0) Y 1/2 W ffi/S (pi 0) Y 1/2 (2.40) 



(no summation implied), which is now interpreted as the Hamiltonian of the dynamical 
system. Explicitly, one finds 

'^ = 0* (24i) 

w aa = w aa . 

Note that this expression coincides with the definition for the Laplacian on a directed graph 
given in [T2|. 

Given the characteristic equation W aj g<£g = A^" 1 one can verify that the ground state is 
now represented by the vector with components 

$) =e -*fl{*)n, (2.42) 

This whole construction can now be interpreted as a quantization procedure by introduc- 
ing a Hilbert space generated by vectors labelled by the states a. Now W is interpreted 
as a Hamiltonian operator, and the evolution equation for cp G as a real time Schrodinger 
equation. 

The analogy to the stochastic quantization for a continuous theory becomes obvious once 
one considers the following points. 

• The Markov property lies at the root of both quantization schemes. 

• Just as Eq. ( 2.32} , also the Fokker-Planck equation ( |2.6} is a master equation. 



The stationary distribution ( 2.34} has the same form as the equilibrium probability in 



Eq. (2.81. 



• The procedure of symmetrizing the matrix W is analogous to the one of bringing the 
FP equation into a Schrodinger-like form leading to the Hamiltonian ( 2.11| . 



Last but not least, 

• the ground states ( |2.12} and \2.42) are the same. 



3 Examples 

This section forms the main body of this article and gives several examples of models which 
are related by the framework of stochastic quantization. We cite some illustrative examples 
from the literature along with examples which we present here for the first time. 
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We start by considering the prototypical example of a zero dimensional field theory 
which becomes supersymmetric quantum mechanics after quantization. We next discuss 
the simple but rich example of the stochastic quantization of a bosonic field. Our treatment 
offers a reinterpretation of the discussion of the quantum Lifshitz model in [2] and allows 
the generalization beyond the free case. We next discuss the relation of the gauged wzw 
model to the strong coupling limit of topologically massive gauge theory which belongs 
to the same class of models even though stochastic quantization cannot be applied directly. 
The resulting limit of topologically massive gauge theory is, unlike the pure Chern-Simons 
theory dynamical. 

Then we move on to the discrete examples of the quantum dimer model and the quantum 
crystal. Here we stress the point that this discrete quantization scheme is the analog of 
stochastic quantization and explore possible interpretations of the results in terms of string 
theory. Finally we sketch the quantization of the six-dimensional Kahler gravity action. 

Recently theories with anisotropic scaling which can be thought of as stemming from a 
stochastic quantization process have appeared in the literature in the context of gauge theory 
membrane actions, and non — Lorentz invariant gravity |13llT4llT5| . 

Another example which we are not discussing here deserves being mentioned, namely 
the connection between d = 3 Chern-Simons theory and d = 4 topological Yang — Mills 
theory via stochastic quantization, as given in [fl6l . 

3.1 From zero dimensions to supersymmetric quantum mechanics 

The simplest example of stochastic quantization is obtained by considering zero dimensional 
quantum field theory. In this case, the field is a map from a point P to the real line, x : P — > R, 
i.e. a variable. The action is a function of this variable, S c i = S c \(x), and the classical partition 
function is given by the integral 



The stochastic quantization is performed by adding a time direction and promoting x to a 
function of time, x : R — > R, t \—> x(t). We can now impose the Langevin equation given in 



where t]{t) is a white Gaussian noise. This system corresponds to a one dimensional Brow- 
nian motion. Using the procedure outlined above, we can write the quantum partition func- 
tion as: 




(3.1) 



Eq. |Z2}: 




(3.2) 



Z = 




dx{t)\ 



(3.3) 
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whence we can directly read the quantum action, which is given precisely by the usual su- 
persymmetric quantum mechanics, 



s q [*,v,fl = / dt \lm 2 +l(s'Mt)}) 2 -mm-ls^[x(t)}mm 



(3.4) 



where tp and xp describe a complex fermion. The corresponding Hamiltonian density is given 
by 

1 „ 1 . . . o 1 .. 

(3-5) 



where n is the conjugate momentum to x(t). Being a supersymmetric theory, it makes sense 
to introduce the charges 



such that the Hamiltonian is obtained as the anticommutator 



n = -{Q,Q} 



(3.6) 
(3.7) 

(3.8) 



The ground state of the theory is annihilated by both charges and can be written as: 

ipo[x] = e - s dM/2. (3.9) 



As expected, the classical partition function factorizes, 

Z d = J^dx e- s *W = Jdx \ip [x}\ 2 = (tp \rpo) L 2 

3.2 Stochastic quantization of a bosonic field 



(3.10) 



As a slightly more complicated example, we treat next the stochastic quantization of a free 
boson. Consider the action for a free boson in d Euclidean dimensions, 



d d x 



difix 1 ) a>(*') 



i = \,2,...,d. 



The Langevin equation, describing the evolution in the fictitious time f is given by 

dtcpihx 1 ) = ^d i d i cp(t > x i ) + ri(t / x i ). 



(3.11) 



(3.12) 



where t](t,x') is a white Gaussian noise, i.e. a stochastic variable whose second momentum 



is the only one that is non-vanishing, see (2.3 1. In this way, f(t, x 1 ) itself becomes a stochastic 
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variable and the expectation value of any functional F[<p] is obtained by averaging over the 
noise: 

<f M>, = \\ -D V F[cp}e-y dtddx ^ 2 , (3.13) 
where the partition function is defined by 

Vt]e - 1 2ldtd"xfj(t,x i )\ 



z 



(3.14) 



It is convenient to change the integration variable from rj to <p. The expression becomes: 



Z 



det 



St] 
Sep 



l]=dt(p(t,x')+Kdjd'cp 



(3.15) 



The Jacobian is easily expressed by introducing two fermionic fields xp(t, x l ) and \p{t, x') such 
that 

5rj 
Sep 



det 



r]=dtf(t,x')+Kdjd'q> 

In this way we can directly read off the (d + 1) -dimensional action: 



(3.16) 



S* +1 [<p,ip,$] = - J dtd d x 



(3.17) 



We can rewrite this result by introducing the superfield O defined by 

0,0) = <p(t,d) +9ip(t,x i ) +xp(t,x i )6 + 6eF(t,x i ) , 



(3.18) 



where 9 and 6 are Grassmann variables and F(t, x 1 ) is an auxiliary bosonic field. One finds 
that 



S d q +1 [f,xp,f] = J d9d6dtd d x [DODO + £ d [<!>]] = J d9d9dtd d x 



DODO + -3,09' O 



(3.19) 

where D and D are the super-covariant derivatives (2.22 1. Note that here, as remarked 
before, C c \ takes the role of a superpotential. In this form the action is manifestly covariant 
under the supersymmetric variations 



S e <p = exp - ipe, 
5 £ ip = e(<p + F), 



S e F = eip + ipe, 
5 £ \p = e(<p-F). 



(3.20) 
(3.21) 



There are different possible approaches to obtain the Hamiltonian description. One pos- 
sibility consists in starting with the Fokker-Planck equation. Alternatively, one can use the 
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supersymmetric structure and define the two supercharges 



Q = f[in (p --d i d t cp) / 



Q = ip(-in cp --d i d' f ) , 



(3.22) 



where EL 



{ Q, Q } = and 



-ij^ is the conjugate momentum to (p. One can easily see that { Q,Q } 



{Q,Q} = 2H, (3.23) 
where Ji is the Hamiltonian operator obtained by a Legendre transformation of the action 



Sq +1 in Eq. (3.17 1. Using the expressions for Q and Q, one can verify that the (bosonic) 
ground state, 

|Y ) = e -!/ d ^^^> |0) , (3.24) 
where |0) is the vacuum annihilated by ip, is in fact a supersymmetric ground state, since 



Q|Y ) = Q|Y ) =0. 



(3.25) 



The bosonic part of the action in Eq. ( |3.17| l has already been conjectured in [17] and studied 
in |2|. The case d = 2 is particularly interesting since the ground state expectation values of 
the d = 3 quantum theory can be evaluated as correlators of a conformal field theory. 

The supersymmetric formulation presented here has the advantage that due to the pres- 
ence of the fermionic part, both Q and Q annihilate the bosonic ground state, which was 
not the case for the quantum Lifshitz model discussed in [2j. In their case, the contribution 
to the Hamiltonian coming from QQ was interpreted as a UV divergent zero-point energy 
and had to be subtracted. Using the stochastic quantization formalism thus allows us to deal 
with more general problems where the boson evolves in a potential and is described by an 
action of the type 



Take for example a massive theory, 



}di<pff<p + V[<p] 



(3.26) 



V[<p] 



m- 



-<p 



(3.27) 



The (d + 1) -dimensional action is found to be 



d d+1 : 



. (3.28) 



This action provides a supersymmetric generalization of the Landau free energy expression 
used to describe the Lifshitz points in the study of liquid crystals ||T8| (a Lifshitz point oc- 
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curs where a high-temperature disordered phase, a spatially uniform ordered phase and a 
spatially modulated phase meet). The supersymmetric charges are found to be 



Q = $ 



-d ; 3> + — <p) , Q = xp 



b K _ ffl 

d,d (P H (P 

<5<p 2 ' r 2 r 



(3.29) 



Another interesting example is provided by starting with a sine-Gordon type action. 



S d d [cp} = f d d x 
In this case, the quantization yields 



^di<p&<p+—CQs(2n<p) 



(3.30) 



St +1 [<p,ip,ip] 



d d+1 : 



\{M) 2 + \ (K9 ! 9> + Asin(27r ( p) x2 



-I/? ( 3t — -yKdid 1 — 7rAcos(27r^) ) i/> , (3.31) 



and the Q and Q operators read: 



O = f[jZ~l 9 ' 9 > " \ sin(27r«p) j , O = f 



~-ldiV<p-^sm{27ap)) . (3.32) 



Since for a precise choice of A and k, the sine-Gordon action corresponds to a massive Dirac 
fermion, we expect this action to be related to the continuum limit of the quantum crystal. It 
was argued in [2] that this action should give rise to a mass gap which is in agreement with 
the results in [3] concerning the discrete model. It is likely that the supersymmetric structure 
we identified here will help in gaining a better understanding of the properties of the system 
away from the Lifshitz point represented by the free boson (quantum dimer). 



3.3 From the gauged wzw model to topologically massive gauge theory 

As the next easiest model to study, let us now consider the gauged WZW model. This ex- 
ample differs a little from the others since we are not using stochastic quantization directly. 
But it belongs to this collection because the two theories (in d, resp. d + 1 dimensions) are 
related by the classical partition function of one being the norm square of the ground state 
wave function of the other. This is the common property of all the examples discussed here. 

We first summarize the WZW model very briefly and refer the reader to the literature 
for details. The basic field g of the wzw model is a map from a 2d Riemann surface E to a 
compact Lie group G. Its basic functional is the action 

M = J^a^'Trig-^g-g-'djg) - iT\g], (3.33) 
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where p is a metric on E and T is the Wess-Zumino term, 



T[§] = l^r j/ 3 ™ m ^g~ X d ig - g- l d jg - g- l d k g, (3.34) 

where B is a 3-manifold such that dB = E. The main property of the WZW model that we will 
make use of in the following is the existence of a holomorphic factorization. This means that 
the partition function Z can be expressed as a square, Z = (f\f), of a holomorphic section / 
of a flat vector bundle over moduli space. 

As it was already noted in Q, a direct stochastic quantization of the WZW action is not 
possible because of the presence of the WZ terrrQ Instead we follow Witten's treatment Ifl9l , 
where the WZW model is gauged in order to derive the existence of the holomorphic factor- 
ization property. In the gauged model, g becomes a section of a bundle X — > E with fiber G 
and structure group Gl x Gr. We take A to be a connection on such a bundle. Let us now 
define 

= M + ^ i d2zTrA ^ ld -S-^ / E d 2 zTrA 2 -A 2 . (3.35) 
Now we can formally define a functional of A, 



Y[A] = J Vge- kI ^> A \ (3.36) 



Y[A] obeys two key equations, 



k 



8A ? 47T 



Ag\ Y[A] = 0, (3.37) 



D ZJTT + ^ D ^z ~ J^F ZZ ) Y[A] = 0, (3.38) 

OAz 4:71 In J 

where the covariant derivative is defined by 

D i u = d i u + [A i ,u]. (3.39) 

These two conditions mean that Y [A] is a holomorphic section and gauge invariant. This can be 
interpreted to mean that Y[A] is a physical state, i.e. a wave function, of 2 + 1 dimensional 
Chern-Simons theory, and thus relates the WZW model to Chern-Simons theory. In [fl9l it 
is shown that the Y[A] we constructed in Eq. ( 3.36| l is indeed the holomorphic section that 



squares to the partition function, 

Z(E) = — f VAW[A}W[A} = ||Y|| 2 . (3.40) 
vol(G) J 



1 The WZ term, being imaginary, would not contribute to the modulus square of the action. As a result, one 
would find at thermal equilibrium, t — > oo, that the system converges to the principal chiral model instead of the 
expected WZW model. 
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We recognize here the by now familiar structure of a theory resulting from stochastic quan- 
tization: the partition function of a classical theory in d dimensions (in this case WZW and 
d = 2) is expressed as the square of the ground state wave function of a quantum theory 
in d + 1 dimensions. All we need to know now is the precise theory whose Hamiltonian 
annihilates Y [A] . The corresponding theory was described in (20l| and is a strong coupling 
limit of topologically massive gauge theory whose action is given by 

S = -jd 3 x Tr[^F 0i F° ! ] + ^J dx T *i AdA ~ l lA "] • ( 3 - 41 ) 

We would now like to point out another similarity with stochastic quantization. Interest- 
ingly enough, the Chern-Simons term is reminiscent of the total derivative term % that 



we dropped in the derivation of the action in Eq. (2.181. When considering a 3-manifold 
Y = E x [0, T], on the one hand the contribution of J Y dxdt^ is J E=a y (S(T) — S(0)), and on 
the other, it is a known fact that on such a manifold, the CS action is equivalent to the (chiral) 
WZW model on the boundary dY |21ll22| . 

Note moreover that in the strong coupling limit we consider, the magnetic component of the 
Yang-Mills term drops out. As a result, the Lagrangian, as it is usually the case for stochastic 
quantization, is not Poincare invariant. Canonically quantizing this model onZxR (where 
E is the Riemann sphere) in the Weyl gauge Aq = 0, one obtains the Hamiltonian^] 

H = J d 2 x Tr[E(x'') 2 ] , i = 1,2, (3.43) 

where 

E fl ' ! '(x) = U a \x) - ^-e^A a Ax) , (3.44) 

OTC ' 

and YI a,l (xi) is the conjugate momentum of Af(x } ): 

[Al{x),Il){y)]=i5 i j5 ab 5{x-y). (3.45) 
In this formalism, we can now understand the two conditions in Eq. ( 3.37} and Eq. (3.38 1: 



• Since Aq appears in the action as a Lagrange multiplier, the gauge choice Aq = is 
implemented by imposing Eq. ( 3.38| as the constraint that the state |Y) be physical. 



• Equation | 3.37| is now understood as Y[A] being the ground state annihilated by the 



2 If instead of taking the strong coupling limit we had considered the complete Yang-Mills term, this would 
have resulted in an extra term in the Hamiltonian, 

H = J d 2 x Tr[£(x') 2 + B(x,) 2 ] , (3.42) 

which does not lead to the ground state we expect. 
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operator E a : 

e " t w = (t^-^) yia1 = - (3 ' 46) 

(Here, we passed to complex coordinates where A z = A\ + 1A2 and A% = A\ — 1A2). 

Being a first order equation, we would like to argue that the last relation can be un- 
derstood as the annihilation of the ground state in a supersymmetric theory. To do so, we 
introduce the following operators: 

Q a (x)=x(x) (v^-^l) =X(x) X E(x), (3.47) 
Q 8 (*)= (~^a+^ A i)x(x) = 2(x) X (x) t (3.48) 
where % is a complex fermion satisfying 

{*(*),x(y)} = *(*-y)- (3-49) 

The anticommutator between Q a and Q a yields 

{Q a (x),Q a (y)} = (V(x)E fl (x) + lLx(x) x ( x )j S(x - y) , (3.50) 

where we used the commutation relation 

[E a (x),E b (y)} = ± T 5 ab 5(x-y). (3.51) 

The anticommutator is a supersymmetric extension of the Hamiltonian of topologically mas- 
sive gauge theory, see Eq. ( 3.43| , which only contains the bosonic part. Using the charges 



Q a (x), Q (x), one can directly derive the supersymmetry transformations for the fields A a 
and X/ which are given by 

5 £ A a = [eQ b + eQ h ,A a }, (3.52) 

S eX = [eQ b + eQ b ,x]- (3-53) 

Going to the Lagrangian formalism, the system is described by the action 

S = -Jd 3 x Tr^F ''] +XX+ J^XX + ^ / dx Tx[AdA - ^iA 3 } . (3.54) 

This can be seen as a strong coupling limit of the supersymmetric topologically massive 
gauge theory studied in 
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3.4 The quantum dimer model 

This example comes from the world of solid state physics. The quantum dimer model is a 
so-called resonating valence bond model and was first introduced by Rokhsar and Kivelson [4] 
as a candidate model for high temperature superconductivity. It starts from the dimer mode^ 
a lattice model living on a bipartite two-dimensional graph. Dimers live on the edges of the 
graph and each vertex can only be touched by one dimer. In the dimer model, the number of 
so-called perfect matchings is counted. The basic move in the dimer model is a plaquette flip, in 
which the dimers of a fully occupied plaquette of the lattice are each turned by one position. 
The quantization of the dimer model is an example of a discrete analog of stochastic quan- 



tization as discussed in Section 2.4 Its Hilbert space is spanned by vectors in one-to-one 



correspondence with the perfect matchings and its quantum Hamiltonian can be expressed 



as 



n = -}(£ !■)(□! + !□)(■! + E + !■)(■!) , (3-55) 



where the black and white squares correspond to the two states in which a fully occupied 
plaquette can be. The sum runs over all flippable plaquettes of a given perfect matching. The 
first two terms are kinetic and flip a plaquette either to the left or the right, the last two terms 
are potential terms and count the number of plaquette moves. The ground state corresponds 
to the equally weighted sum over all perfect matchings. 

In the case of the dimer model on the hexagonal lattice, there exists a one-to-one map 
between perfect matchings and the configurations of an idealized three-dimensional crystal 
corner. The Hamiltonian in Eq. ( 3.55} can be generalized to describe a (Markov) growth 



process as detailed in the next section. 



3.5 Quantum crystal melting 



Let us now discuss the example of quantizing the melting crystal corner [3|, which is inter- 
esting from the point of view of string theory. This is as well an example of a discrete analog 



of stochastic quantization as discussed in Section 2.4 



The melting crystal corner is often mapped to the problem of stacking cubes in an empty 
corner of 3D space. The growth rules are the following. A cube can be added to a configu- 
ration if three of its sides will touch either the wall or other cubes. This leads to a minimum 
energy configuration without any free-standing cubes. The partition function of the melting 
crystal corner takes the following form: 



Z 



1 



E « =Ylj^aW' (356) 

3d partitions n=l \ 1 J 



3 For the detailed definitions see for example |24| 
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where the rightmost expression is the so-called MacMahon function. 

Crystal melting is an extension of the dimer model discussed in the last section in the 
following sense. Via a rhombus tiling, the configurations of the dimer model on the hexag- 
onal lattice are in one-to-one correspondence with the configurations of the melting crystal 
corner. While in the dimer model, only the overall number of configurations are counted, 
the model of the melting crystal corner contains more information, it also keeps track of the 
number of cubes of each configuration. The quantum Hamiltonian of the melting crystal can 
be expressed as 



H = 



-7(x:i")(n| + P)("l + Ev^PX°l + 4l")("l) , (3-57) 

which acts on the Hilbert space generated by the orthonormal basis of the configurations 
| a). The sum runs over all places in a given configuration where a cube can be added or 
removed. The ket ■) represents a place where a cube can be removed, while (□) represents 
a place where a cube can be added. Once one takes the limit q — > 1 in ( 3.57} , one recovers 
the Hamiltonian of the quantum dimer model. The ground state is in fact the unique zero 
energy ground state, fulfilling 

H (ground) = 0, (3.58) 

and has the form 

(ground) =Y j q N{a)/1 \oi) . (3.59) 

a 

Note the normalization for the wave function, 

(ground|ground) = y^ #boxes = Z . (3.60) 

a 

Once we identify 

q = e-8 s , (3.61) 

it turns out that the closed string partition function of the A-model with the Calabi-Yau 
X = C 3 as target space corresponds exactly to the partition function of the melting crys- 



tal corner (3.561 USUI. F° r the crystal, the above identification (3.611 translates into a lat- 
tice spacing of g s . The crystal melting configurations can be roughly mapped to (multiple) 
blow-ups X of the original Calabi-Yau X. Moreover, in [25J it was shown that a crystalline 
structure similar to the one corresponding to C 3 which we used, can be constructed for any 
toric Calabi-Yau manifold, based on the corresponding quiver gauge theory. The partition 
function of this crystalline model will then correspond to the partition function of the topo- 
logical A-model on this toric Calabi-Yau modulo possible wall crossings (see e.g. J26J). This 
means that also our quantization scheme can be generalized to any toric Calabi-Yau. All that 
is needed is the state graph whose nodes are crystal configurations. Its arrows mark which 
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configurations are related by adding or removing an atom. The quantization procedure does 
not depend on the details of the crystal. Consider the state graph for the melting of a crystal 
made of A different atoms. To each configuration oc we associate a weight 



w(ct) 



(3.62) 



i=l 



where < q, < 1, i = 1, . . . , A are parameters and tij(a) is the number of atoms of species 



i in the configuration a. Then the Laplacian W in Eq. (2.41 1 becomes the matrix with entries 



a and /3 are neighbours, 



E 



, V w(a) 

fa V v ' 



if a = p, 



(3.63) 



7 neighbour of t 

otherwise. 



One can verify that the vector 



(ground) = E EN/ 

a \ i=l 



n,-(«)/2 



\X) 



(3.64) 



is a zero energy ground state for W and its norm square reproduces the classical crystal 
melting partition function 

(ground | ground) = }w{ct) = Z ~ Z top (modulo wall crossings) . (3.65) 



Let us now consider possible interpretations of the quantized crystal in terms of string 
theory. By adding a time evolution to the statistical system, we created a quantum theory 
with one dimension more than the system we started from. 

From the string theory point of view, this means that instead of the six dimensional topo- 
logical A-model, we should be looking at a seven dimensional theory. An immediate candi- 
date would be topological M-theory or an effective version thereof. Given that the new time 
dimension is not treated on the same footing as the original six dimensions, we cannot be 
looking at M-theory itself, but rather at a non-Lorentz invariant limit. Unfortunately, topo- 
logical M-theory has not yet been clearly defined. Its status is best discussed in [8], where 
its connection to Hitchin's functionals is described. 

Given that the classical configurations of the melting crystal are related to geometries, 
it would be interesting to try to find a completely geometric interpretation of the Hamilto- 
nian ( 3.57| , which is the Laplacian on the graph of states (geometries). Its kinetic term creates 
"neighboring" geometries from a given configuration. It would be interesting to see if also 
the potential term has a precise geometric meaning. 
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Concerning the interpretation of the quantum crystal, another observation can be made 
that we consider to be one of the central conceptual points of this paper. If we identify the 
three-dimensional crystal configurations with multiple blow-ups of C 3 , we find that the 
Hilbert space generated by the configurations has a natural interpretation as mini-superspace 
and the ground state can be seen as a wavefunction of the Universe. 

Looking at formula ( 3.581 1, 

H\Y Q )=0, (3.66) 

one is strongly reminded of the Wheeler-De Witt equation. The WDW equation also has the 
form Ti-c^ = 0, where Hq is the Hamiltonian associated to the Einstein-Hilbert action, and 
Y is the wave function of the universe. The WDW equation is in fact nothing else but a zero 
energy Schrodinger equation and the crucial information needed is the initial quantum state 
which is then evolved. The wave function Y is defined on superspace, an infinite dimensional 
space of all possible geometries and matter configurations. In practice, one usually works in 
the finite dimensional mini-superspace, in which all but a few degrees of freedom are frozen 
out. 

The sum over classical crystal melting configurations corresponds on the A-model side 
to a sum over "quantum Kahler geometries", or to be more precise, a sum over ideal sheaves, 
which are torsion-free sheaves with vanishing first Chern class [6|. In this framework, the 
mini-superspace we are looking at is thus the moduli space of all ideal sheaves J\A shea f over 
the original manifold X. 

In our case however, it is not obvious how to find the natural classical gravity theory 
whose quantization leads to Ji (i.e. the analog of the Einstein-Hilbert action the usual WDW 
equation descends from). Let us therefore concentrate on the meaning of Ji. As we have 
stressed already, H is the Laplacian on the space of configurations. The condition in Eq. ( 3.66} 



can be understood as requiring Yo to be a harmonic function in mini-superspace. This is 
strongly reminiscent of the results in [27] concerning quantum gravity in 2 + 1 dimensions. 
In this case, the Einstein-Hilbert gravity turns into a Hamiltonian system on Teichmiiller 
space after the ADM reduction. In particular, for genus g = 1, the WDW equation takes the 
form 

v / A|Y)=0, (3.67) 

where A is the Laplacian on the torus moduli space. Two observations can be made about 
this: 

• (quantum) gravity in 2 + 1 dimensions is special because there are no propagating 
degrees of freedom. It is therefore not surprising that this case is very close to our 
construction which is related to topological strings. 



• The square root in ( |3.67 1 is reminiscent of the fact that in supersymmetric systems the 



ground state is found by imposing Q Yq = 0. The operator Q is in some sense close 
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to being the "square root" of the Hamiltonian, since U = \{Q,Q}. Also the discrete 
quantum crystal can be supersymmetrized (by enlarging the space of states) and the 
ground state then satisfies a first order difference equation as opposed to the second 
order equation ( 3.66| . Apart from this, Eq. ( 3.67\ results in the same ground state as 



Eq. ( pT66) . 



There is yet another way of looking at the classical crystal configurations, as adopted 
in (25| , namely as BPS bound states of D-branes (DO and D2 branes bound to a single D6 
brane). Removing an atom of the crystal corresponds here to adding a specific combination 
of DO and D2 brane charges. In this picture, the quantum Hamiltonian ( 3.57\ jumps thus 
between configurations with different D-brane charges. 

In the above paragraphs, we have outlined possible directions to pursue in order to ar- 
rive at a meaningful interpretation of the quantum crystal on the string theory side of the 
correspondence. We leave a more thorough investigation of these ideas for future research. 

3.6 Kahler gravity? 

In the previous section we proposed a possible interpretation of our quantum crystal con- 
struction in terms of a Wheeler-De Witt equation. To make this connection more precise 
we would need an appropriate theory of gravity to quantize. A possible candidate is a 6- 
dimensional form theory of gravity, namely Kahler gravity, which describes variations of 
the Kahler structure on a complex manifold M |7] . Also this example is interesting from the 
point of view of string theory, since it provides the target space description of the topological 
A-model. Following the logic in |5l |6], its stochastic quantization should provide a theory 
that corresponds to (or is at least in the same universality class and thus shares important 
properties with) the quantum crystal. A relation of this theory to topological M-theory is 
therefore to be expected. Here we sketch a possible approach leading to the stochastic quan- 
tization of this theory. 

The action of Kahler gravity is given by 

SKaehler [K] = J ^K-^dK + AKAK, (3.68) 

where K is a variation of the (complexified) Kahler form on M, and d c = d — d. The Kahler 
gravity action is invariant under gauge transformations of the form 

6 a K = da — d cf (K A a) , (3.69) 

where a is a 1-form on M, such that d cf a. = 0. The equations of motion for Kahler gravity 
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theory have the form 



dK + d' 



(KAK) = 0. 



(3.70) 



We can decompose K into massless and massive modes, 



K = x + 




x G H U (M,C), 



(3.71) 



where x G H 1 ' 1 (M / C) represents the Kahler moduli, which are not integrated over, and 
7 G Q 3 (M) contains the massive modes of K. Using the above decomposition we can write 
the Kahler gravity action without non-local terms as follows: 



Lagrangian D-branes of the A-model are charged under 7, implying that these branes are 
sources for K and hence modify the integral of K on 2-cycles which link them. This also 
implies that the partition function of the A-model depends non-perturbatively on the choice 
of a cohomology class in H 3 (M) as well as on x G H 2 (M). 

Let us now implement the stochastic quantization procedure in the supersymmetric ap- 
proach. The seven-dimensional theory can be written as 



where D and D are the super-covariant derivatives, * is the Hodge star on the 6d manifold, 
and T is the superfield defined by 



where tp, ip, F are all three-forms in space. Note that x contains the Kahler moduli and is not 
promoted to a superfield. Just as we have seen above, the classical action now takes the role 
of a superpotential. 

The above derivation is purely formal: the properties of the seven dimensional action 
have to be understood, and it is not obvious that the stochastic quantization of Kahler gravity 
is the most direct way towards a description of the seven dimensional topological M-theory. 
On the other hand we believe that some of the general features that we encountered in the 
previous examples (quantum theory in d + 1 dimensions whose ground state reproduces 
the partition function of the classical theory in d dimensions, supersymmetry in the extra 
dimension) are to be expected for such a theory. 




(3.72) 




(3.73) 



r = 7 + 9tp + ip6 + 66F, 



(3.74) 
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4 Conclusions 

In this article we have discussed a rather widespread, even though often unrecognized, 
scheme to relate a classical to a quantum field theory. In condensed matter physics, this 
method of quantization gave rise to the (discrete) quantum dimer model and its general- 
ization, the quantum crystal. In the continuous case, this approach goes under the name of 
stochastic quantization. In topological field theories it has appeared for example in relating 
the WZW model to the Chern-Simons action. Its main features are the following. 

1 . The classical theory lives on a d-dimensional manifold E, while the (d + 1 ) -dimension- 
al (quantum) theory lives on the manifold Y = E x [0, T]. 

2. In the Hamiltonian description, the quantum theory admits a zero-energy ground state 
wavefunction Yo. Its norm square reproduces the classical partition function Z c \\ 

(Y |Y ) = Z cl . (4.1) 

3. The quantum theory can be interpreted as a Markov process that at equilibrium con- 
verges to the minima of the classical action. In this sense, the quantum theory describes 
the same physics in the T — > oo limit as the classical theory. 

4. There is a natural supersymmetric structure in the extra dimension. 

Based on point number 2 above, one might even venture to say that every theory that 
admits a holomorphic factorization of its partition function might have a (d + 1) -dimensional 
partner whose ground state square results in the partition function of the original theory. 

In this note, we have supplied a number of examples of models which are as such widely 
known, while their relation to each other by stochastic quantization on the other hand is 
largely unrecognized. 

We have done two things: 

1. Having recognized a common principle in seemingly unrelated models, we are able 
to overcome some of the problems encountered in previous works, notably by making 
use of the inherent supersymmetry. 

2. Since this quantization scheme underlies the quantum crystal, which is deeply related 
to topological strings, we would like to propose it as an approach towards topological 
M-theory. 

We have sketched some ideas for interpreting the quantum crystal on the string theory 
side of the correspondence, exploring interpretations in the WDW framework and in terms 
of a non-Lorentz invariant limit of topological M-theory. Much work remains to be done in 
both directions and a more thorough study of these ideas is left for future research. 
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